Plasma-Derived Exosome Proteins as Novel Diagnostic and Prognostic Biomarkers in Neuroblastoma Patients

Neuroblastoma (NB) is the most common extracranial solid tumor during infancy, causing up to 10% of mortality in children; thus, identifying novel early and accurate diagnostic and prognostic biomarkers is mandatory. NB-derived exosomes carry proteins (Exo-prots) reflecting the status of the tumor cell of origin. The purpose of this study was to characterize, for the first time, the Exo-prots specifically expressed in NB patients associated with tumor phenotype and disease stage. We isolated exosomes from plasma specimens of 24 HR-NB patients and 24 low-risk (LR-NB) patients at diagnosis and of 24 age-matched healthy controls (CTRL). Exo-prot expression was measured by liquid chromatography–mass spectrometry. The data are available via ProteomeXchange (PXD042422). The NB patients had a different Exo-prot expression profile compared to the CTRL. The deregulated Exo-prots in the NB specimens acted mainly in the tumor-associated pathways. The HR-NB patients showed a different Exo-prot expression profile compared to the LR-NB patients, with the modulation of proteins involved in cell migration, proliferation and metastasis. NCAM, NCL, LUM and VASP demonstrated a diagnostic value in discriminating the NB patients from the CTRL; meanwhile, MYH9, FN1, CALR, AKAP12 and LTBP1 were able to differentiate between the HR-NB and LR-NB patients with high accuracy. Therefore, Exo-prots contribute to NB tumor development and to the aggressive metastatic NB phenotype.


Introduction
Neuroblastoma (NB) is the most common extracranial solid tumor during infancy, with a median age at diagnosis of 17 months.Although classified as a rare disease, it represents a worldwide emergency, causing up to 10% of mortality in children.NB shows notable heterogeneity with regard to histology and clinical behavior, ranging from low-risk (LR-NB) localized tumors to high-risk (HR-NB) disease, characterized by an aggressive metastatic phenotype, resistance to treatment, and fatal relapse occurrence.Risk stratification demands the highest accuracy, as it determines the therapeutic treatment.For NB patients, risk assessment is based on age at diagnosis, stage of the disease, chromosomal aberrations Cells 2023, 12, 2516 2 of 28 and amplification of MYCN, a transcriptional factor considered the major oncogenic driver.However, the current therapeutic stratification, based on such clinical and molecular risk factors, does not allow to discriminate patients with similar clinical-pathological parameters who receive the same treatment, despite showing markedly different clinical courses.The limits of the existing classification system could be overcome by a deeper investigation of the biology of the tumor.Thus, the challenge remains the identification of additional and more accurate prognostic markers to improve risk stratification and direct a personalized treatment.
Despite aggressive therapy, the five-year survival rate of HR-NB patients remains poor.Moreover, HR-NB tumors are often unresectable as they infiltrate surrounding areas, hindering tissue biopsies.Nowadays, the molecular characterization of neoplastic diseases allows to overcome the limitations of the cell morphology-based classification of solid tumors, improving cancer diagnosis and treatment [1].The acquired experience in molecular profiling highlights some limitations, especially in terms of sampling methods [2].Although they are still the main source of biological material for diagnostic purposes, tissue biopsies entail high costs, high risk of adverse effects because of the invasive nature of the technique and procedural complications resulting in inadequate material retrieval and lead to increased false-negative or false-positive results that affect treatment decisions [3,4].Moreover, tissue sampling does not capture all of the spatial and temporal heterogeneity of tumors, preventing the identification of aggressive and therapy-resistant cellular subclones that are responsible for disease progression and recurrence [5].
Recently, interest has grown in liquid biopsies, which represent a surrogate of the tumor and can provide complementary or even additional information for tissue biopsies.The investigation of biological fluids is based on a minimally invasive, less expensive and easily accessible method that can be applied for cancer diagnosis and prognosis, treatment selection, disease monitoring and assessment of relapse occurrence.Liquid biopsies are a source of circulating tumor cells, circulating tumor DNA, exosomes, and nano-sized vesicles (30-150 nm) derived from the endocytic pathway and released in the extracellular environment after the fusion of multivesicular bodies with the cell plasma membrane [6].The release of exosomes is regulated by specific stimuli, and their cargo reflects the producing cell [7].Indeed, exosomes are defined as the bioprint of the cell of origin and, thus, the study of their content allows to gain accurate information of primary tumor features [8][9][10].Indeed, exosomes represent a source for both diagnostic and prognostic markers in several human diseases [8][9][10][11][12][13][14][15].To gain insight into the characteristics of NB for improving diagnosis and treatment, we performed proteomic analysis of the exosomes released in the blood of NB patients.In particular, in order to identify molecules with potential value for improving risk stratification, we compared the exosomal protein (Exo-prot) profiles of LR-NB and HR-NB patients at the onset.Our results identified the Exo-prots and biological pathways involved in tumor development and progression and associated with the acquisition of aggressive tumoral traits, integrating the data of our previous study demonstrating that exosomal microRNAs (exo-miRNAs) are associated with HR-NB patients' response to treatment [8].

Study Cohort and Blood Sample Collection and Processing
The study cohort included Italian NB patients classified into different risk groups: HR-NB (n = 24) and LR-NB (n = 24).Age-matched blood samples of healthy donors (n = 24) were used as the control (CTRL).The main clinical features of the study cohort are reported in Table 1.The procedures for patient enrollment were carried out according to the proper guidelines and in adherence with the ethical principles of the Declaration of Helsinki.Blood samples were provided by the BIT-Gaslini Biobank, which centralizes the NB specimens derived from Italian AIEOP (Associazione Italiana Emato-Oncologia Pediatrica) centers.Written informed consent was provided from the parents or legal guardian of each patient considered for the study.Blood samples were collected in EDTA tubes and processed within Cells 2023, 12, 2516 3 of 28 24 h.Processing included centrifugation at 1200× g for 10 min at room temperature (RT) to collect plasma.Plasma was stored at −80 • C or used immediately for exosome isolation.

Exosome Isolation and Proteomic Analysis
Exosome isolation was performed from 500 µL of plasma utilizing an exoRNeasy Serum/Plasma Midi Kit (Qiagen Italia, Milan, Italy) with a modified version of the protocol that we have previously described [8], which enables to isolate vesicles with diameter and surface markers specific to exosomes.Protocol modifications were added starting from the washing step with XWP buffer, performed twice by centrifuging samples at 500× g and 5000× g for 5 min at RT, respectively.The isolated exosomes were lysed, reduced and alkylated directly on the filter by adding to 100 µL of boiling lysis buffer 6 M Guanidine, 10 mM TCEP, and 4 mM CAA in 100 mM Tris pH 8 and heating the columns in a thermomixer at 100 • C for 10 min and 1000 rpm.The lysed samples were centrifuged at 4000× g for 5 min and left to digest for 2 h with the addition of 0.6 µg of LysC (FUJIFILM Wako Pure Chemical Corporation, Milano, Italy).The filters were then washed with 200 µL of Tris 25 mM, pH 8, by centrifuging at 4000× g for 5 min.The total sample was further digested by adding 0.6 µg of Trypsin (Promega, Milano, Italy), followed by mixing and incubating at 37 • C overnight.The digested samples were loaded onto StageTips [16].The resulting peptides were completely dried in a SpeedVac centrifuge at 30 • C, resuspended in 2% ACN and 0.1% formic acid and analyzed by a nano-UHPLC-MS/MS system using an Ultimate 3000 RSLC coupled to an Orbitrap Fusion Tribrid mass spectrometer (Thermo Scientific, Milano, Italy).

Proteomic Set-Up
Elution was performed with an EASY spray column (75 µm × 50 cm, 2 µm particle size, Thermo Scientific, Milano, Italy) at a flow rate of 250 nL/min with a 100 min non-linear gradient of 7-45% solution B (80% ACN, 20% H2O, 5% DMSO and 0.1% FA).MS analysis was performed in DDA mode.Orbitrap detection was used for MS1 measurements at a resolving power of 120 K in a range between 375 and 1500 m/z, with an AGC target of 400,000 and a maximum injection time of 50 ms.MS/MS spectra were acquired in the linear ion trap (rapid scan mode) after collision-induced dissociation (CID) fragmentation at a collision energy of 35% and an AGC target of 2000 for up to 300 ms.A 1.5 s cycle time was performed for data-dependent MS/MS analysis, during which precursors with a charge range of 2-4 were selected for activation in order of abundance.Quadrupole isolation with a 1.8 m/z isolation window was used, and dynamic exclusion was enabled for 25 s.
Raw data were processed with MaxQuant [17] software, version 1.6.4.0.A false discovery rate (FDR) of 0.01 was set for the identification of proteins, peptides and PSM (peptide-spectrum match).For peptide identification, a minimum length of 6 amino acids was required.The Andromeda engine, incorporated into MaxQuant software, was used to search the MS/MS spectra against the Uniprot human database and additional human database (release UP000005640_9606 April 2018).In the processing, acetylation (Protein N-Term), oxidation (M) and deamidation (NQ) were selected as variable modifications, and the fixed modification was carbamidomethylation (C).The quantification intensities were calculated by the default fast MaxLFQ algorithm with the activated option "match between runs."

Bioinformatic Analysis
Differential expression analysis was carried out using the SMAGEXP tool suite [18].Exo-prots with a Benjamini-Hochberg p-value lower than 0.05 and a log fold change >0.58 or <−0.58 were considered statistically significantly deregulated.Exo-prots with ≤30% missing values were retained for the analysis to reduce the bias introduced by imputation.The random forest (RF) and quantile regression imputation of left-censored data (QRILC) methods were used for missing value imputation using the METimp web tool [19].The RF method uses an iterative imputation approach that trains an RF on observed values and then predicts the missing values [20].The QRILC method performs the imputation of left-censored missing data, using random draws from a truncated distribution with parameters estimated by quantile regression [19].Protein-protein interaction (PPI) network and pathway analyses were carried out on selected statistically significantly deregulated Exo-prots using STRING-DB software (version 12.0) [21].The PPI enrichment p-value was used to assess whether a group of proteins have more interactions among themselves than what would be expected from a random set of proteins of a similar size, drawn from the genome.We considered significant a PPI enrichment p-value lower than 0.05.Pathway analysis was used to identify the biological processes mainly regulated by differentially expressed Exo-prots.An FDR lower than 0.05 identified significantly enriched ontology terms and pathways.Venn diagrams were plotted using the InteractiVenn web-based tool [22].

Statistical Analysis
A two-proportion Z-test based on the detectable and missing values between the NB and CTRL or HR-NB and LR-NB groups was used to determine whether the two proportions were different from each other.A Z-score <−2.6 or >2.6 were considered to demonstrate statistically significant differences.Age group and sex were considered potential confounders, and differences between the NB patients and CTRL were tested using Fisher's exact test.A Fisher p-value lower than 0.05 was considered statistically significant.Receiver operating characteristic (ROC) curves were plotted to visually display the discriminating power of each significantly modulated Exo-prot across groups.Exo-prot expression was used to generate ROC curves.The area under the ROC curve (AUC) was computed to quantitatively assess the discriminating performance of each Exo-prot.ROC curves and statistical calculations were performed using GraphPad Prism 8.0 software (GraphPad Software, La Jolla, CA, USA).
The combined discriminatory power of a selection of statistically significantly represented Exo-prots across the NB and CTRL groups was assessed by generalized linear models using the GLMNET R package [23].The group was set up as the response variable, Exo-prot expression as the explanatory variable, the elastic net mixing parameter at 1.0, and the minimum mean square error as the gamma value.The best gamma value was assessed using the leave-one-out cross-validation technique.

ROC Data Validation by Gene Expression
Diagnostic or prognostic Exo-prots identified by ROC analysis were validated by measuring the expression of the corresponding genes in the primary tumor tissue at the onset.Tissue biopsies from 10 LR and 10 HR-NB patients and four non-oncologic control tissues were homogenized with the Tissue Lyser II instrument and RNA was extracted with an RNeasy Mini Kit (Qiagen Italia, Milan, Italy).The RNA quality was assessed with a Nano RNA Assay on the Agilent 2100 Bioanalyzer (Agilent Technologies Spa, Milan, Italy).RNA quantification was performed with the QIAxpert System (Qiagen Italia, Milan, Italy).RNA was reverse-transcribed with SuperScript™ II Reverse Transcriptase (Thermo Fisher Scientific, Milan, Italy) and the expression levels of the FN1, AKAP12, MYH9, CALR, NCL, NCAM, VASP, LGALS3BP, DCN and LUM genes were assessed by RTqPCR (real-time quantitative PCR) with specific TaqMan assays (Thermo Fisher Scientific, Milano, Italy).Each sample was repeated in triplicate on a MicroAmp Fast Optical 96-well reaction plate run on the ViiA7 Real-Time PCR System (Thermo Fisher Scientific, Milan, Italy).GAPDH was used as the housekeeping gene for data normalization.The relative expression values between sample groups were calculated using the 2 (−DCT) method.Statistical analysis and plots were carried out using GraphPad Prism 8.0 for Windows.

Protein Cargo Profiling of the NB and CTRL Subjects
To identify the potential biomarkers and key biological processes associated with tumor development in NB patients, we performed proteomic analysis of the exosomes isolated from the plasma samples of the NB patients and age-matched CTRL provided by the BIT-Gaslini biobank.Table 1 summarizes the clinical and molecular characteristics of the patients and CTRL.A flowchart summarizing the main steps of the sample collection and bioinformatic analyses is shown in Figure 1.Specifically, plasma samples were collected from 24 LR-NB patients (stage 1, 2, 4 s and MYCN non-amplified tumors), 24 HR-NB patients (age >18 months and stage 4 tumor or MYCN amplified) and 24 CTRL children.The association between sex (male vs. female) or age group (<18 months vs. >18 months) and the subject group (CTRL) or NB-HR vs. LR-NB was analyzed to assess the potential confounding effects of these characteristics.Fisher's exact test did not show a significant association between sex or age group and the subject group (p > 0.05), thereby excluding the possibility that sex and age group might be confounding factors in this study.
Plasma-derived exosomes from each sample were isolated, and protein cargo was profiled by high-resolution mass spectrometry coupled with liquid chromatography.The analytical data generation identified 458 exosomal proteins (Exo-prots) detectable in at least one NB or CTRL subject.

Analysis of Detectable Exo-Prots Based on the Distribution of Missing Values in the NB and CTRL Subjects
Since missing values can adversely affect statistical analysis [24], the distribution of the Exo-prot missing values between the NB patients and CTRL was investigated to assess whether detectable Exo-prots might provide new insights into disease development.The Z-test analysis of the detectable value proportions showed that 38 Exo-prots were exclusively expressed in the NB patients, whereas 83 were exclusively expressed in the CTRL subjects (z-score >2.6 or <−2.6;Supplementary Table S1).A heatmap based on the expression of these 121 Exo-prots visually shows a clear association with the NB patients and CTRL subjects, visually confirming the results of the Z-test analysis of proportions (Figure 2A).Among the exclusively detectable Exo-prots in the NB patients, we identified the MYC target nucleophosmin (NPM1) and the actin interacting protein zyxin (ZYX), known for their role in tumor proliferation and invasiveness [25,26].

Analysis of Detectable Exo-Prots Based on the Distribution of Missing Values in the NB an CTRL Subjects
Since missing values can adversely affect statistical analysis [24], the distribution the Exo-prot missing values between the NB patients and CTRL was investigated to a sess whether detectable Exo-prots might provide new insights into disease developmen The Z-test analysis of the detectable value proportions showed that 38 Exo-prots we exclusively expressed in the NB patients, whereas 83 were exclusively expressed in th CTRL subjects (z-score >2.6 or <−2.6;Supplementary Table S1).A heatmap based on th expression of these 121 Exo-prots visually shows a clear association with the NB patien and CTRL subjects, visually confirming the results of the Z-test analysis of proportio (Figure 2A).Among the exclusively detectable Exo-prots in the NB patients, we identifie the MYC target nucleophosmin (NPM1) and the actin interacting protein zyxin (ZYX known for their role in tumor proliferation and invasiveness [25,26]. Network analysis was carried out on the 121 Exo-prots to help interpret the resu and assess the biological function of these proteins (Figure 2B).The network analys displayed significantly more interactions than expected for a random set of proteins of similar size drawn from the genome (protein-protein interaction (PPI) enrichme p-value <0.05).These findings establish a connection between significantly modulate Exo-prots and indicate that these proteins are biologically connected as a group.Pathwa analysis is a well-known bioinformatic tool that is used to explore the biological pr cesses and pathways associated with a list of differentially expressed genes/proteins u ing curated ontologies [27].Network analysis was carried out on the 121 Exo-prots to help interpret the results and assess the biological function of these proteins (Figure 2B).The network analysis displayed significantly more interactions than expected for a random set of proteins of a similar size drawn from the genome (protein-protein interaction (PPI) enrichment p-value <0.05).These findings establish a connection between significantly modulated Exo-prots and indicate that these proteins are biologically connected as a group.Pathway analysis is a well-known bioinformatic tool that is used to explore the biological processes and pathways associated with a list of differentially expressed genes/proteins using curated ontologies [27].Pathway analysis identified 141 significantly enriched biological processes a pathways in the NB patients (Supplementary Table S2).A selection of the most sign cantly enriched terms are listed in Table 2.These results show that the identifi Exo-prots in the NB patients might act as cancer-associated pathways, such as stress sponse, immune system regulation, inflammation, cell motility and cytoskeletal re rangements.Pathway analysis identified 141 significantly enriched biological processes and pathways in the NB patients (Supplementary Table S2).A selection of the most significantly enriched terms are listed in Table 2.These results show that the identified Exo-prots in the NB patients might act as cancer-associated pathways, such as stress response, immune system regulation, inflammation, cell motility and cytoskeletal rearrangements.Missing values cannot be directly handled by bioinformatics tools.Imputation is a bioinformatic technique used to substitute missing values with defined observations.We applied RF or QRILC, two well-known methods [19], for imputing missing values to subsequently identify deregulated Exo-prots.

Differentially Expressed Exo-Prots in the NB Patients
Exo-prots showing at most 30% missing values were retained for imputation, while the rest were filtered out (See Section 2.4).The two imputation methods successfully substituted missing values with numeric values, thus defining two distinct datasets, hereafter referred to as the RF and QRILC datasets.Principal component analysis based on the protein expression profile of 24 HR-NB, 24 LR-NB and 24 CTRL subjects was performed in order to assess the presence of outlier samples.The results indicate that subject CTRL_18 was a potential outlier exclusively in the QRLIC dataset, but was not in the RF dataset.Therefore, we concluded that QRLIC imputation could have had a substantial impact on the proteomic profile of this specific subject.In order to identify specific deregulated Exo-prot biomarkers in the NB patients, we compared the Exo-prot expression profile of the NB and CTRL subjects using the SMAEXP tool in the RF and QRILC datasets.In the RF dataset, differential expression analysis identified 38 upregulated and 90 downregulated Exo-prots in the NB patients compared to the CTRL subjects (Supplementary Table S3), whereas in the QRILC dataset, 32 upregulated and 92 downregulated Exo-prots were identified (Supplementary Table S4).A Venn diagram of the statistically significantly modulated Exo-prots in the two datasets showed that 17 upregulated and 64 downregulated Exo-prots in the NB patients were shared by the two datasets.Meanwhile, 21 upregulated and 26 downregulated Exo-prots were exclusive to the RF dataset and 15 upregulated and 28 downregulated Exo-prots were exclusive to the QRILC dataset, respectively (Figure 3A).A complete list of the 17 upregulated and 64 downregulated Exo-prots in the NB patients identified in both datasets is reported in Table 3. Heatmap visualization based on the expression values of the 81 Exo-prots commonly regulated in the RF and QRILC datasets shows a clear separation between the NB and CTRL groups (Figure 3B,C).Notably, among the upregulated Exo-prots, we identified proteins that positively regulate tumor development and progression, such as neural cell adhesion molecule (NCAM) [28], nucleolin (NCL) [29] and galectin-3-binding protein (LGALS3BP) [30].Among the downregulated Exo-prots, we identified proteins with oncosuppressive functions, such as lumican (LUM) [31], decorin (DCN) [32] and the vasodilator stimulated phosphoprotein (VASP) [33].Network analysis was carried out on the 17 upregulated and 64 downregulated Exo-prots (Figure 3D), showing their strong interactions (PPI enrichment p-value <0.05).These findings establish various interactions between significantly modulated Exo-prots in Table 2 and indicate that these proteins are biologically connected as a group.Pathway analysis identified 153 significantly enriched biological processes and pathways in the NB patients (Supplementary Table S5).A selection of the most significantly enriched terms are listed in Table 4.These results show that the NB patients were characterized by aberrant expression of proteins involved mainly in inflammation and complement activation, immune response, cell proliferation and apoptosis and extracellular matrix (ECM) interactions.biologically connected as a group.Pathway analysis identified 153 significantly enrich biological processes and pathways in the NB patients (Supplementary Table S5).A sel tion of the most significantly enriched terms are listed in Table 4.These results show th the NB patients were characterized by aberrant expression of proteins involved mainly inflammation and complement activation, immune response, cell proliferation a apoptosis and extracellular matrix (ECM) interactions.

Analysis of the Detectable Exo-Prots in HR-NB or LR-NB Based on the Distribution of Missing Values
To identify detectable Exo-prots in the HR-NB patients, we compared the Exo-prot expression profiles of the HR-NB and LR-NB patients.A Z-test was used to identify Exoprots whose distribution of missing values was significantly different between the HR-NB and LR-NB patients.This analysis showed that 63 Exo-prots were exclusively expressed in the HR-NB patients, whereas 52 were exclusively expressed in the LR-NB patients (Z-score >2.6 or <−2.6;Supplementary Table S6).A heatmap based on the expression of these 115 Exo-prots showed a clear discrimination between the outcome groups, confirming the results of the Z-test analysis of proportions (Figure 2C).Among the exclusively detectable Exo-prots, we identified well-known inflammation biomarkers such as CPR, SAA1/SAA2 and PTX3 [34,35].Network analysis showed statistically significant interactions among 115 Exo-prots (PPI enrichment p-value <0.05; Figure 2D).Pathway analysis identified 205 significantly enriched processes in the NB patients (Supplementary Table S7).A selection of the most significantly enriched terms are listed in Table 5.These results show that the identified Exo-prots in the HR-NB patients are mainly involved in immune responses, inflammation, cell migration, apoptosis and cytoskeletal organization.

Differentially Expressed Exo-Prots in the HR-NB Patients
To identify deregulated Exo-prots in the HR-NB patients, we analyzed the Exo-prot expression profiles between the HR-NB and LR-NB patients using the SMAEXP tool in the RF or QRILC dataset.In the RF dataset, the analysis identified 98 upregulated and 70 downregulated Exo-prots in the HR-NB patients compared to the LR-NB patients (Supplementary Table S8), whereas in the QRILC dataset, 58 upregulated and 59 downregulated Exo-prots were identified (Supplementary Table S9).A Venn diagram of the deregulated Exo-prots in the two datasets showed that 46 upregulated and 46 downregulated Exo-prots in the HR-NB patients were common between the two datasets.On the contrary, 52 upregulated and 24 downregulated Exo-prots were exclusive to the RF dataset and 12 upregulated and 13 downregulated Exo-prots were exclusive to the QRILC dataset, respectively (Figure 4A).A complete list of 46 upregulated and 46 downregulated Exo-prots in the HR-NB patients identified in both datasets is reported in Table 6.Heatmap visualization of the expression values shows a clear separation between the HR-NB and LR-NB groups (Figure 4B,C).In the HR-NB patients, we observed the upregulation of proteins that positively regulate cell migration and metastasis: Myosin-9 (MYH9) [36,37], fibronectin (FN1) [38] and latenttransforming growth factor-beta-binding protein 1 (LTBP1) [39].Among the downregulated Exo-prots, we identified proteins exerting antitumoral effects: Calreticulin (CALR) [40] and A-kinase anchor protein 12 (AKAP12) [41][42][43].Network analysis was carried out on the 46 upregulated and 46 downregulated Exo-prots in the HR-NB patients (Figure 4D).These findings establish statistically significant interactions among the deregulated Exo-prots in Table 6 and indicate that these proteins are biologically connected as a group (PPI enrichment p-value <0.05).Pathway analysis identified 482 significantly enriched biological processes and pathways in the HR-NB patients (Supplementary Table S10).A selection of the most significantly enriched terms is reported in Table 7.These results show that the Exo-prots modulated in the HR-NB patients are involved not only in inflammation, immune response and cell proliferation, but also specifically in cell migration, adhesion and motility, as well as in angiogenic processes.Our data suggest that Exo-prots may promote the metastatic properties of HR-NB tumors.In order to assess the sensitivity and specificity of the potential biomarkers, we performed ROC analysis on the selected Exo-prots among those found to be significantly differentially expressed in the NB and HR-NB patients.Analysis was performed using both the RF and QRILC datasets.NCAM, NCL, LGALS3BP, LUM, DCN and VASP were studied for their potential diagnostic value in NB, whereas MYH9, FN1, LTBP1, CALR and AKAP12 were investigated for their prognostic value in HR-NB.The area under the curve (AUC) was used to measure the degree of separability of the Exo-prots between the NB and CTRL or HR-NB and LR-NB patients.The results show that NCAM (RF: AUC 0.83; QRILC: AUC 0.83, both p < 0.0001; Figure 5A), NCL (RF: AUC 0.85, QRILC: AUC 0.83, both p < 0.0001; Figure 5B), LUM (RF: AUC 0.77, p = 0.0002/QRILC: AUC 0.77, p = 0.0002; Figure 5C) and VASP (RF: AUC 0.8, QRILC: AUC 0.85, both p < 0.0001; Figure 5D) had the highest AUCs and p-values, indicating a robust and significant diagnostic value for NB patients.Meanwhile, LGALS3BP (RF: AUC 0.65, p = 0.03/QRILC: AUC 0.65, p = 0.03; Figure 5E) and DCN (RF: AUC 0.75, p = 0.0005/QRILC AUC 0.71, p = 0.003; Figure 5F) have a low but significant AUC, indicating a lower diagnostic value.We also investigated whether the combination of Exo-prots could improve the discriminatory power of these markers taken singularly.We applied the generalized linear model (GLMNET) algorithm (see Material and Methods).Our findings show that the discriminatory power of the combined model was higher than that obtained by each biomarker taken singularly in both datasets (RF: AUC 0.97, QRILC: AUC 0.97, both p < 0.0001; Figure 5G).
Comparing the HR-NB and LR-NB samples, ROC analysis evidenced that MYH9 (RF: AUC 0.95, QRILC: AUC 0.92, both p < 0.0001; Figure 6A), FN1 (RF: AUC 0.89, QRILC: AUC 0.89, both p < 0.0001; Figure 6B) and CALR (RF: AUC 0.94, QRILC: AUC 0.94, both p < 0.0001; Figure 6C) had high AUCs and low p-values, indicating a robust and significant prognostic value.LTBP1 (RF: AUC 0.78, p = 0.0007/QRILC AUC 0.71, p = 0.01; Figure 6D) and AKAP12 (RF: AUC 0.74; p = 0.005/QRILC: AUC 0.84, p < 0.0001; Figure 6E) had lower AUCs, but were still significant, indicating a lower prognostic value.Similarly, as performed between the NB and CTRL groups, we assessed whether the combination of Exo-prots could improve the discriminatory power of these markers taken singularly.We applied the GLMNET algorithm and showed that, again, the discriminatory power of the combined model was higher than that obtained by each single marker in both datasets (RF: AUC 0.99, QRILC: AUC 0.97, both p < 0.0001; Figure 6F).Comparing the HR-NB and LR-NB samples, ROC analysis evidenced that MYH9 (RF: AUC 0.95, QRILC: AUC 0.92, both p < 0.0001; Figure 6A), FN1 (RF: AUC 0.89, QRILC: AUC 0.89, both p < 0.0001; Figure 6B) and CALR (RF: AUC 0.94, QRILC: AUC 0.94, both p < 0.0001; Figure 6C) had high AUCs and low p-values, indicating a robust and significant prognostic value.LTBP1 (RF: AUC 0.78, p = 0.0007/QRILC AUC 0.71, p = 0.01; Figure 6D) and AKAP12 (RF: AUC 0.74; p = 0.005/QRILC: AUC 0.84, p < 0.0001; Figure 6E) had lower Exo-prots could improve the discriminatory power of these markers taken singularly.We applied the GLMNET algorithm and showed that, again, the discriminatory power of the combined model was higher than that obtained by each single marker in both datasets (RF: AUC 0.99, QRILC: AUC 0.97, both p < 0.0001; Figure 6F).In conclusion, these results indicate that the expression of exosomal NCAM, NCL, LUM and VASP has evident diagnostic value in discriminating NB patients from CTRL subjects, both taken singularly or in combination, suggesting that these non-invasive Exo-prot biomarkers may be used as novel diagnostic tools and provide new insights in the mechanisms of neuroblastoma pathogenesis.Meanwhile, the expression of MYH9, FN1, LTBP1, CALR and AKAP12 is effective in discriminating HR-NB from LR-NB patients, indicating that these non-invasive biomarkers may be used to detect patients at risk of developing aggressive tumor phenotypes.In conclusion, these results indicate that the expression of exosomal NCAM, NCL, LUM and VASP has evident diagnostic value in discriminating NB patients from CTRL subjects, both taken singularly or in combination, suggesting that these non-invasive Exoprot biomarkers may be used as novel diagnostic tools and provide new insights in the mechanisms of neuroblastoma pathogenesis.Meanwhile, the expression of MYH9, FN1, LTBP1, CALR and AKAP12 is effective in discriminating HR-NB from LR-NB patients, indicating that these non-invasive biomarkers may be used to detect patients at risk of developing aggressive tumor phenotypes.

Validation of the Key Selected Proteins by Gene Expression
To validate the diagnostic and prognostic Exo-prots identified by ROC analysis, we measured the expression of the corresponding coding genes in tissue biopsies of the NB primary tumors of 10 HR and 10 LR cases.We used tissue samples of age-matched nononcologic subjects as the control.NCL and NCAM were confirmed as diagnostic markers, significantly overexpressed in NB tumors compared to control samples (p < 0.01 and p < 0.05, respectively).The diagnostic power of DCN and LUM was confirmed as well, with a significant downregulation of both markers in NB patients (p < 0.001 and p < 0.01, respectively) (Figure 7A).VASP did not show any significant modulation between the two biological groups considered, while LGALS3BP, contrary to the exosome data, was downregulated in the NB primary tumor tissues (Figure 7A).From a prognostic point of view, significant FN1 overexpression (p < 0.05) was confirmed in the HR-NB compared to LR-NB patients, while MYH9 and LTBP1 did not report any significant modulation (Figure 7B).In agreement with the proteomic data, AKAP12 and CALR were significantly decreased in HR disease (p < 0.05 and p < 0.01, respectively) (Figure 7B).
primary tumors of 10 HR and 10 LR cases.We used tissue samples of age-matched non-oncologic subjects as the control.NCL and NCAM were confirmed as diagnostic markers, significantly overexpressed in NB tumors compared to control samples (p < 0.01 and p < 0.05, respectively).The diagnostic power of DCN and LUM was confirmed as well, with a significant downregulation of both markers in NB patients (p < 0.001 and p < 0.01, respectively) (Figure 7A).VASP did not show any significant modulation between the two biological groups considered, while LGALS3BP, contrary to the exosome data, was downregulated in the NB primary tumor tissues (Figure 7A).From a prognostic point of view, significant FN1 overexpression (p < 0.05) was confirmed in the HR-NB compared to LR-NB patients, while MYH9 and LTBP1 did not report any significant modulation (Figure 7B).In agreement with the proteomic data, AKAP12 and CALR were significantly decreased in HR disease (p < 0.05 and p < 0.01, respectively) (Figure 7B).

Discussion
In the present study, we characterized the protein cargo of the exosomes derived from patients with LR-NB and HR-NB tumors and age-matched control subjects to identify novel non-invasive diagnostic or prognostic biomarkers.
We applied LC-MS-based proteomics, a technique that has increased the knowledge about the protein content of extracellular vesicles in recent years [44].In vitro proteomics studies on neuroblastoma-derived exosomes have been conducted to determine their role in tumor progression and metastasis formation [45][46][47][48].However, these studies have not been carried out on patient-derived samples, thus limiting their translational applicability to the clinical settings.
The proteomic profile of each patient and CTRL subject was analyzed with two distinct bioinformatic approaches-the distribution of missing values and differential expression analysis.Distribution analysis is expected to detect Exo-prots in which missing values are exclusively detectable in one of the biological groups, whereas differential expression analysis aims to assess the deregulation of Exo-prot expression in a large majority of samples.We considered that missing value distribution analysis could provide useful insights into tumor development because Exo-prots that are undetectable or detectable only in NB patients could actually have diagnostic value and potentially a biological meaning.Distribution analysis can identify distinct Exo-prots from those found via differential expression analysis because no filter is applied to the number of missing values.On the contrary, the differential expression analysis was based on Exo-prots with a maximum of 30% missing data that underwent missing value imputation.Thus, the results of both analyses were complementary.
Analysis of the distribution of missing values was instrumental to identify Exo-prots in NB patients involved in the stress response, immune system regulation, inflammation, cell motility and cytoskeletal rearrangements.Among the exclusively detectable Exo-prots in the NB patients, we identified the nucleolar protein NPM1 as being overexpressed in hematologic and solid tumors [49] and ZYX involved in focal adhesion and in the regulation of actin cytoskeleton with a role in cell motility and tumor invasion [50].
The comparison between HR-NB and LR-NB patients identified the C-reactive protein (CRP), recognized as a biomarker of inflammation, whose elevated serum level has been shown to be associated with poor prognosis in many cancers and in NB [51]; serum amyloid A1/2 (SAA1/SAA2), which have been found to be prognostic in many solid tumors, and whose expression, in combination with high level of CRP, can stratify patients with high-risk early-stage melanoma [34]; pentraxin 3 (PTX3), which is reported to be associated with the PI3K/AKT/mTOR signaling pathway to induce tumor cell proliferation and apoptosis and whose overexpression results in poor survival outcomes in many tumors [52].Although the missing value distribution analysis suggested the pro-tumoral activity of deregulated Exo-prots, missing data can negatively affect statistical power [24].Thus, imputation of the missing values is usually performed to obtain robust results [24].Since RF and QRILC are indicated among the most accurate imputation methods [19], we decided to use both approaches to perform missing value imputation.
Differential expression analysis was instrumental in identifying Exo-prots in the NB patients involved in tumor-promoting processes, such as inflammation and complement activation, immune response, cell proliferation and apoptosis and ECM interaction.Among the upregulated Exo-prots, we identified NCAM, an adhesion molecule known to be highly expressed in NB tumors and associated with increased metastatic capacity [28], and LGALS3BP, a glycosylated protein that has migratory-promoting activity [30].Moreover, NCL was upregulated in the NB patients; this phosphoprotein regulates several biological processes, including cell proliferation and apoptosis, DNA and RNA metabolism, ribosomal RNA maturation and ribosome assembly [29].Of note, a high expression of NCL has recently been demonstrated to be an unfavorable independent prognostic factor in NB [53].Hence, the demonstration of NCL overexpression in exosomes released in the circulation of patients might be relevant to better understand its role in driving NB tumor progression.
Among the downregulated proteins, we identified two proteoglycans belonging to the small leucine-rich proteoglycan (SLRP) family: LUM and DCN.The former exerts an oncosuppressive function by impairing metalloproteinase 14 (MMP-14)-mediated cell migration [31], while the latter can inhibit tumor cell growth and migration and, simultaneously, induce the autophagic process of stromal cells, thus exerting a double action that results in a sharp oncosuppressive function [32].We also reported the downregulation of VASP, an assembly regulatory protein interacting with ZYX, and regulating cell adhesion and migration [33].Interestingly, the absence of VASP/ZYX interaction impairs cell adhesion and promotes cell migration and invasion.Moreover, VASP affects the phosphorylation of focal adhesion kinase (FAK), which has anti-apoptotic properties, thus being involved in the regulation of the apoptotic process [33].Tumor-derived exosomes carry on their surface specific tumor markers.Disialoganglioside GD2 is a marker associated with neuroectoderm-derived cancers and has been identified on the surface of exosomes isolated from melanoma [54] and neuroblastoma patients [8].Among the Exo-prot cargo, we were not able to observe specific enzymes involved in the ganglioside metabolic pathway differentially expressed between the NB and CTRL subjects.This result was also confirmed by pathway analysis, which did not reveal the enrichment of ganglioside synthesis/catabolism processes.
The comparison between HR-NB and LR-NB patients identified Exo-prots upregulated and downregulated in patients with HR-NB disease.Network and pathway analyses revealed that such proteins are connected and act in common biological systems, exerting a major function in pathways that are not only involved in oncogenesis, such as inflammation, cell proliferation and immune response regulation, but also in the acquisition of tumor aggressive traits, such as loss of cell adhesion, migration and angiogenesis.In particular, we observed the upregulation of MYH9, which regulates cell polarity and cytoskeleton [36,37] and has been shown to enhance cell migration in breast cancer cells [37].Indeed, the exosome-mediated transport of MYH9 is able to directly promote macrophage migration, resulting in increased macrophage infiltration at the tumor site and, consequently, in cancer cell metastasis formation [36,37].The significantly deregulated Exo-prots in the HR-NB patients included molecules involved in ECM interaction, a key aspect for the establishment of cell-to-cell communication, with a major impact on the definition of the tumor microenvironment required for tumor growth and dissemination.Interestingly, among the upregulated Exo-prots in the HR-NB patients, we identified FN1 and LTBP1.The former is a glycoprotein that exerts its oncogenic role mainly by upregulating cell proliferation and inhibiting apoptosis [38,55] and by regulating angiogenesis, contributing to the direction of blood vessel formation and stimulating cell differentiation toward an endothelial phenotype [38].LTBP1 is an extracellular matrix protein that enhances cell migration and invasion and promotes the epithelial-to-mesenchymal transition (EMT), a process at the basis of the acquisition of cell motility and invasion properties, leading to metastasis formation [39].In esophageal squamous cell carcinoma, a positive correlation between FN1 and LTBP1 has been demonstrated [39].Indeed, while LTBP1 was mainly expressed in tumor tissues, FN1 was highly expressed in fibroblasts of the stroma, suggesting a potential correlation between the overexpression of the two genes, which synergistically act to induce EMT.Among the downregulated proteins, we identified CALR, a protein localized in the endoplasmic reticulum and surveilling proper protein folding, cell adhesion, integrin signaling and antigen presentation [40].CALR exposed on the surface of tumor cells can trigger phagocytosis and initiate the anticancer immune response [40].Hence, the downregulation of CALR can have oncogenic effects because of the impairment of physiological homeostasis and the reduced stimulation of immune surveillance against cancer cells [40].We also observed downregulation of AKAP12, a scaffold protein belonging to the A-kinase anchoring protein (AKAP) family and regulating signal transduction cascade.In particular, AKAP members drive the compartmentalization of cAMP-dependent protein kinases (PKAs) to different subcellular locations to ensure the specificity of PKA-mediated signal transduction [56].AKAP12 exerts a tumor-suppressive function through the negative modulation of angio-genesis, cell proliferation and migration by preventing PKC activation [57].Thus, reduced expression of AKAP12 may result in the activation of pro-tumoral and metastatic pathways.
Importantly, ROC analysis demonstrated that the exosomal levels of NCAM, NCL, LUM and VASP have significant diagnostic power in discriminating NB patients from CTRL subjects.High levels of NCAM and NCL have already been described as unfavorable prognostic factors in NB [53,58], and herein we reported, for the first time, their significant upregulation within the exosomes derived from plasma samples of NB patients at diagnosis.In addition, our results indicated that MYH9, FN1, CALR, AKAP12 and, with a lower performance, LTBP1 expression can discriminate HR-NB and LR-NB patients with high specificity and sensitivity.Interestingly, FN1 overexpression has been associated with mesenchymal phenotype [59], which confers chemoresistance to NB cells.Thus, the exosome-mediated delivery of such proteins that we reported to be overexpressed in HR-NB patients could be involved in therapy-resistant mechanisms.Our results showed that the combination of all Exo-prot diagnostic and prognostic markers outperformed each single marker.Therefore, measuring the whole panel of such molecules will provide higher specificity and sensitivity to the analysis.Importantly, the overexpression of NCL and NCAM and the downregulation of DCN and LUM in the NB patient-derived exosomes compared to the CTRL samples was confirmed also in primary tumor tissues.Consistency between the proteomics data and tumor gene expression was also reported for the upregulation of FN1 and the downregulation of CALR and AKAP12 in the HR-NB patients compared to the LR-NB cases.The lack of correspondence between exosomes and tumor tissues for VASP, LGALS3BP, MYH9 and LTBP1 could be due to the fact that these vesicles are loaded through active mechanisms [60].Thus, the differential expression specifically observed in the exosomes may be due to the cargo selection, which does not necessarily reflect the tissue gene expression levels.Functional experiments and a prospective clinical trial should be performed to increase our knowledge of the biological mechanisms involved in NB formation and development and to robustly assess the applicability of Exo-prots for improving NB patient diagnosis and prognosis.

Conclusions
We characterized the Exo-prot cargo of HR, LR-NB and CTRL subjects and identified significantly differentially expressed molecules.We assessed the diagnostic significance of specific deregulated Exo-prots and we identified highly accurate markers differentiating NB vs. CTRL samples and HR vs. LR-NB patients.The present study is the first, to the best of our knowledge, to characterize the protein content of exosomes derived from the plasma samples of NB patients.The results showed that these vesicles contain differentially expressed Exo-prots that are mainly involved in cancer-associated pathways and in the acquisition of aggressive tumor features, thus playing a major role in driving oncogenesis and tumor progression.Importantly, we demonstrated that a few specific Exo-prots are able to discriminate, with high accuracy, NB patients from CTRL samples, whereas others are able to precisely distinguish HR-NB from LR-NB patients.These results provide evidence of the potential relevance of liquid biopsies to integrate the diagnosis of NB tumors and improving risk stratification to help refine patient treatment.

Supplementary Materials:
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cells12212516/s1.Table S1: List of the Exo-prots exclusively expressed in the NB patients compared to the CTRL samples; Table S2: Significantly enriched processes associated with the Exo-prots exclusively detectable in the NB and CTRL subjects; Table S3: Differentially expressed Exo-prots between the NB and CTRL subjects in the RF dataset; Table S4: Differentially expressed Exo-prots between the NB and CTRL subjects in the QRILC dataset; Table S5: Significantly enriched processes associated with the Exo-prots regulated in the comparison between the NB and CTRL samples; Table S6: List of Exo-prots exclusively detectable in the HR-NB or LR-NB samples; Table S7: Significantly enriched processes associated with the Exo-prots exclusively detectable in the HR-NB or LR-NB samples; Table S8: Differentially expressed Exo-prots between the HR-NB and LR-NB subjects in the RF dataset; Table S9: Differentially expressed Exo-prots between

Figure 1 .
Figure 1.Flowchart of the analysis.Flow chart of the NB (n = 48) and CTR (n = 24) plasma sample analysis.NB included LR-NB (n = 24) and HR-NB (n = 24) patients.The diagram reports the different steps of the bioinformatic analysis applied to the proteomic data, which includes the distribution analysis of missing data, data filtering based on the number of missing values, the application of two different imputation methods (RF and QRILC), differential expression analyses, the selection of the Exo-prots commonly differentially expressed between the RF and QRILC datasets and pathway and network analyses.

Figure 2 .
Figure 2. Distribution analysis of the detectable and missing expression values across samples.(A) Heatmap showing the significantly differentially represented Exo-prots (n = 121) between the NB and CTRL samples according to the Z-test statistics based on the distribution of missing values.Exo-prots with Z-test values larger than 2.6 or lower than −2.6 were considered to have a significantly different distribution of missing values between groups.The heatmap was drawn using global coloring.Detectable values are colored in red and missing values in grey.A dendrogram is reported above the plot.(B) Protein-protein interaction network among the 121 differentially represented Exo-prots.A network was built using the STRING-DB tool.The line thickness indicates the strength of data support.A medium confidence of 0.4 was set up as the minimum required interaction score.(C) Heatmap showing the significantly differentially represented Exo-prots (n = 115) between the HR-NB and LR-NB samples according to the Z-test statistics.Exo-prots with Z-test values larger than 2.6 or lower than −2.6 were considered to have a significantly different distribution of missing values across groups.The heatmap was drawn using global coloring.Detectable values are colored in red and missing values in grey.A dendrogram is reported above the plot.(D) Protein-protein interaction network among the 115 differentially represented Exo-prots.A network was built using String-DB.The line thickness indicates the strength of the data support.A medium confidence of 0.4 was set up as the minimum required interaction score.

Figure 3 .
Figure 3. Significantly differentially expressed Exo-prots between the NB and CTRL samples Venn diagram showing the number of commonly and exclusively up-and downregu Exo-prots in the RF and QRILC datasets.(B,C) Heatmaps showing the expression values o deregulated Exo-prots identified in the RF (panel B) and QRILC (panel C) datasets.The NB CTRL group labels are shown above the plot.A color key is reported in the top left part of the (D) Protein-protein interaction network showing the interactions among the 81 differentially resented Exo-prots between the NB and CTRL samples and shared by the RF and QRILC data A network was built using the STRINGg-DB tool.The Exo-prots are shown as nodes.The width indicates the strength of data support.A medium confidence of 0.4 was set up as the mum required interaction score.

Figure 3 .
Figure 3. Significantly differentially expressed Exo-prots between the NB and CTRL samples.(A) Venn diagram showing the number of commonly and exclusively up-and downregulated Exo-prots in the RF and QRILC datasets.(B,C) Heatmaps showing the expression values of the deregulated Exo-prots identified in the RF (panel B) and QRILC (panel C) datasets.The NB and CTRL group labels are shown above the plot.A color key is reported in the top left part of the plot.(D) Protein-protein interaction network showing the interactions among the 81 differentially represented Exo-prots between the NB and CTRL samples and shared by the RF and QRILC datasets.A network was built using the STRINGg-DB tool.The Exo-prots are shown as nodes.The line width indicates the strength of data support.A medium confidence of 0.4 was set up as the minimum required interaction score.

Figure 4 .
Figure 4. Significantly differentially expressed Exo-prots between the HR-NB and LR-NB samples.(A) Venn diagram showing the number of commonly and exclusively up-and downregulated Exo-prots in the RF and QRILC datasets by comparing the HR-NB and LR-NB samples.(B,C) Heatmaps showing the expression values of the deregulated Exo-prots identified in the RF (panel B) or QRILC (panel C) dataset.The HR-NB and LR-NB group labels are shown above the plot.A color key is reported in the top left part of the plot.(D) Protein-protein interaction network showing the interactions among the 92 differentially represented Exo-prots.Network analysis showing the interactions among the 92 commonly differentially expressed proteins between the HR-NB and LR-NB samples and shared by the RF and QRILC datasets.The network was built using the STRING-DB tool.Exo-prots are shown as nodes.The line width indicates the strength of data support.A medium confidence of 0.4 was set up as the minimum required interaction score.

Figure 4 .
Figure 4. Significantly differentially expressed Exo-prots between the HR-NB and LR-NB samples.(A) Venn diagram showing the number of commonly and exclusively up-and downregulated Exoprots in the RF and QRILC datasets by comparing the HR-NB and LR-NB samples.(B,C) Heatmaps showing the expression values of the deregulated Exo-prots identified in the RF (panel B) or QRILC (panel C) dataset.The HR-NB and LR-NB group labels are shown above the plot.A color key is reported in the top left part of the plot.(D) Protein-protein interaction network showing the interactions among the 92 differentially represented Exo-prots.Network analysis showing the interactions among the 92 commonly differentially expressed proteins between the HR-NB and LR-NB samples and shared by the RF and QRILC datasets.The network was built using the STRING-DB tool.Exo-prots are shown as nodes.The line width indicates the strength of data support.A medium confidence of 0.4 was set up as the minimum required interaction score.

3. 6 .
ROC Analysis of the Diagnostic Value of the Exo-Prots in the NB and CTRL Subjects and the Prognostic Significance of the Exo-Prots in the HR-NB and LR-NB Patients

Cells 2023 , 29 Figure 5 .
Figure 5. Results of the ROC analysis to assess the diagnostic value of the selected differentially expressed Exo-prots between the NB and CTRL subjects.ROC analysis was performed to assess the diagnostic value of the differentially expressed Exo-prots between the NB and CTRL samples.Analysis was carried out using the expression values from the RF or QRILC datasets.ROC curves for NCAM (A), NCL (B), LUM (C), VASP (D), LGALS3BP (E) and DCN (F) singularly and the combination of all markers (G) are reported.Sensitivity (%) is shown on the Y-axis, 100%; specificity (%) is shown on the X-axis.Areas under the curves (AUCs), coincidence intervals (CIs) and p-values are reported for each graph.The red dotted line corresponds to AUC = 0.5, representing a random classifier.p-Values lower than 0.05 are considered statistically significant.Red and blue labels refer to the upregulated and downregulated Exo-prots in the NB vs. CTRL samples, respectively.The gene symbol and dataset name are reported above each plot.

Figure 5 .
Figure 5. Results of the ROC analysis to assess the diagnostic value of the selected differentially expressed Exo-prots between the NB and CTRL subjects.ROC analysis was performed to assess the diagnostic value of the differentially expressed Exo-prots between the NB and CTRL samples.Analysis was carried out using the expression values from the RF or QRILC datasets.ROC curves for NCAM (A), NCL (B), LUM (C), VASP (D), LGALS3BP (E) and DCN (F) singularly and the combination of all markers (G) are reported.Sensitivity (%) is shown on the Y-axis, 100%; specificity (%) is shown on the X-axis.Areas under the curves (AUCs), coincidence intervals (CIs) and p-values are reported for each graph.The red dotted line corresponds to AUC = 0.5, representing a random classifier.p-Values lower than 0.05 are considered statistically significant.Red and blue labels refer to the upregulated and downregulated Exo-prots in the NB vs. CTRL samples, respectively.The gene symbol and dataset name are reported above each plot.

Figure 6 .
Figure 6.Results of the ROC analysis to assess the prognostic value of the selected differentially expressed Exo-prots between the HR-NB and LR-NB samples.ROC analysis was performed to assess the prognostic value of selected Exo-prot in discriminating between the HR-NB and LR-NB patients.Analysis was carried out using the expression values from the RF and QRILC datasets.ROC curves for MYH9 (A), FN1 (B), CALR (C), LTBP1 (D) and AKAP12 (E) singularly and the combination of all markers (F) are reported.Sensitivity (%) is shown on the Y-axis, 100%; specificity (%) is shown on the X-axis.Areas under the curves (AUCs), coincidence intervals (CIs) and p-values are reported for each graph.The red dotted line corresponds to AUC = 0.5, representing a random classifier.p-Values lower than 0.05 are considered statistically significant.Red and blue labels refer to upregulated and downregulated Exo-prots in the HR-NB vs. LR-NB patients, respectively.The gene symbol and dataset name are reported above each plot.

Figure 6 .
Figure 6.Results of the ROC analysis to assess the prognostic value of the selected differentially expressed Exo-prots between the HR-NB and LR-NB samples.ROC analysis was performed to assess the prognostic value of selected Exo-prot in discriminating between the HR-NB and LR-NB patients.Analysis was carried out using the expression values from the RF and QRILC datasets.ROC curves for MYH9 (A), FN1 (B), CALR (C), LTBP1 (D) and AKAP12 (E) singularly and the combination of all markers (F) are reported.Sensitivity (%) is shown on the Y-axis, 100%; specificity (%) is shown on the X-axis.Areas under the curves (AUCs), coincidence intervals (CIs) and p-values are reported for each graph.The red dotted line corresponds to AUC = 0.5, representing a random classifier.p-Values lower than 0.05 are considered statistically significant.Red and blue labels refer to upregulated and downregulated Exo-prots in the HR-NB vs. LR-NB patients, respectively.The gene symbol and dataset name are reported above each plot.

Figure 7 .LFigure 7 .
Figure 7. RTqPCR validation of the genes encoding the diagnostic or prognostic Exo-prots identified by ROC analysis.Transcript levels were compared between the HR-NB (n = 10), LR-NB (n = 10) tissues and CTRL subjects (A) and between the HR-NB and LR-NB patients (B).The graph shows the mean of 2 (−DCT) values, against the mean of the DCt of the reference gene GAPDH.The results are shown as a box plot and are expressed as the mean normalized gene expression values.The boxes show the values that fall between the 25th and 75th percentiles, the horizontal lines represent

Table 1 .
Clinical features of the NB and CTRL subjects included in this study.

Table 2 .
Selection of the significant pathways of deregulated Exo-prots in the NB patients according to analysis of the distribution of missing values.
GO BP, KEGG and Reactome pathway analyses were carried out with STRING-DB on Exo-prots exclusively expressed in the NB or CTRL samples.The pathways are listed by increasing false discovery rate values for each ontology.a Name of the ontology defining a pathway.GO BP stands for Gene Ontology biological process.KEGG stands for Kyoto Encyclopedia of Genes and Genomes.b Official name of a GO biological process, KEGG or Reactome pathway.c Number of identified genes.d FDR (false discovery rate) shows p-values corrected for multiple testing within each category using the Benjamini-Hochberg procedure.Values lower than 0.05 are considered significant.

Table 3 .
Relative expression of the commonly regulated Exo-Prots between the NB and CTRL samples in the RF and QRILC datasets.
3.30 × 10 −10Exo-prot expression profile was evaluated in the NB and CTRL samples, and comparative analysis of the expression differences between the two groups was carried out by a two-sample Student's t-test after imputation by the RF or QRILC methods.Exo-prots are listed by decreasing FC value.a Name of the protein-coding genes.b Fold change values expressed as log2 after the RF imputation method.FC values greater than 0.67 are reported.c p-Value adjusted for FDR (false discovery rate).Adjusted p-values lower than 0.05 are considered significant.d Fold change values expressed as log2 after the RF imputation method.FC values greater than 0.67 are reported.e p-Value adjusted for FDR (false discovery rate).Adjusted p-values lower than 0.05 are considered significant.

Table 4 .
Selection of significant pathways of the Exo-prots regulated in the comparison between the NB and CTRL samples.GO BP, KEGG and Reactome pathway analyses were carried out with STRING-DB on the Exo-prots differentially expressed in the NB vs. CTRL samples.Pathways are listed by increasing false discovery rate values for each ontology.
a Name of the ontology defining a pathway.GO BP stands for Gene Ontology biological process.KEGG stands for Kyoto Encyclopedia of Genes and Genomes.b Official name of a GO biological process, KEGG or Reactome pathway.c Number of identified genes.d FDR (false discovery rate) shows p-values corrected for multiple testing within each category using the Benjamini-Hochberg procedure.Values lower than 0.05 are considered significant.

Table 3 .
Relative expression of the commonly regulated Exo-Prots between the NB and C samples in the RF and QRILC datasets.

Table 5 .
Selection of significant pathways of the Exo-prots exclusively detectable in the HR-NB or LR-NB patients.
GO, KEGG and Reactome pathway analyses were carried out with STRING-DB on Exo-prots exclusively expressed in the HR-NB or LR-NB samples.Pathways are listed by increasing false discovery rate values for each ontology.a Name of the ontology defining a pathway.GO BP stands for Gene Ontology biological process.KEGG stands for Kyoto Encyclopedia of Genes and Genomes.b Official name of a GO biological process, KEGG or Reactome pathway.c Number of identified genes.d FDR (false discovery rate) shows p-values corrected for multiple testing within each category using the Benjamini-Hochberg procedure.Values lower than 0.05 are considered significant.

Table 6 .
Relative expression of the commonly regulated Exo-prots between the HR-NB and LR-NB patients in the RF and QRILC datasets.

Table 6 .
Relative expression of the commonly regulated Exo-prots between the HR-NB and LR-NB patients in the RF and QRILC datasets.

Table 6 .
Cont.Exo-prot expression profiles were evaluated in the HR-NB and LR-NB samples, and comparative analysis of the expression differences between the two sample groups was carried out by a two-sample Student's t-test after imputation by the RF or QRILC method.Exo-prots are listed by decreasing FC value.a Name of the protein-coding genes.b Fold change values expressed as log2 after the RF imputation method.FC values greater than 0.67 are reported.c p-Value adjusted for FDR (false discovery rate).Adjusted p-values lower than 0.05 are considered significant.d Fold change values expressed as log2 after the QRILC imputation method.FC values greater than 0.67 are reported.e p-Value adjusted for FDR (false discovery rate).Adjusted p-values lower than 0.05 are considered significant.

Table 7 .
Selection of significant pathways of the Exo-prots regulated in the comparison between the HR-NB and LR-NB samples.

Table 7 .
Cont.GO BP, KEGG and Reactome pathway analyses were carried out with STRING-DB on the Exo-prots differentially expressed in the HR-NB and LR-NB samples.Pathways are listed by increasing false discovery rate values for each ontology.a Name of the ontology defining a pathway.GO BP stands for Gene Ontology biological process.KEGG stands for Kyoto Encyclopedia of Genes and Genomes.b Official name of a GO biological process, KEGG or Reactome pathway.c Number of identified genes.d FDR (false discovery rate) shows p-values corrected for multiple testing within each category using the Benjamini-Hochberg procedure.Values lower than 0.05 are considered significant.